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Abstract: Two-Higgs-doublet models (2HDM) are simple extensions of the Standard 
Model (SM) where the scalar sector is enlarged by adding a weak doublet. As a result, the 
Higgs potential depends in general on several free parameters which have to be carefully 
chosen to give predictions consistent with the current precision data. We consider a 2HDM 
invariant under a twisted custodial symmetry and depending only on three extra parameters 
beyond the SM ones. This model can naturally features an inverted mass spectrum with a 
light pseudoscalar state and a heavy SM-like Higgs boson. We thoroughly analyze direct 
and indirect constraints and present a few unconventional though promising signatures at 
the LHC. 
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1. Introduction 

Even in the absence of any definitive experimental evidence, the mechanism of spontaneous 
symmetry breaking |2|, |3| remains the most promising candidate to consistently explain 
the existence of massive weak gauge bosons. Following Occam's principle, this mechanism 
is implemented within the Standard Model (SM) using a single SU{2)l Higgs doublet. 
Such a minimal solution displays three important advantages. First, it guarantees that 
a zero mass for the photon is not a mere accident Second, it implies that the scalar 
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potential is invariant under an SU{2)l x S'C/(2)/j global symmetry, spontaneously broken 
down to a "custodial" SU (2) ||5| . This custodial symmetry naturally ensures that the 
one-loop quantum corrections to the successful tree-level relation 



where gL and gy are the coupling constant associated with the gauge groups SU{2)l and 
U{1)y, vary at most logarithmically with the Higgs boson mass. Finally, it naturally 
explains the absence of tree-level Flavour Changing Neutral Currents (FCNC) [^. 

Nevertheless, more complex possibilities involving additional and/or larger representa- 
tions may a priori be considered. Motivations for such enlarged Higgs sectors fall in general 
into two (not mutually exclusive) categories. Either they can be associated with extended 
symmetries, such as Supersymmetry and Grand Unification Theory (GUT) groups, or 
they can be justified by phenomenological arguments, such as the observed predominance 
of matter over anti-matter in our Universe. 

The generic two-Higgs-doublet model, where the SM scalar sector is enlarged with only 
one additional SU{2)l doublet, appears as an attractive laboratory. Though relatively 
simple, it can display a rich phenomenology. Specific realisations of this model can also be 
predicted, or required, by Beyond the Standard Model (BSM) theories such as the Minimal 
Supersymmetric SM (MSSM) or S'O(10)-invariant GUTs. However, contrary to the SM 
scalar sector, the most general 2HDM Higgs potential does not guarantee = 0, neither is 
invariant under custodial symmetry nor under CP symmetry, and leads to large FCNC. In 
order to naturally explain the absence of large corrections to precisely measured quantities, 
additional symmetries in the Higgs sector of the theory must therefore be imposed. In 
particular, the existence of a custodial symmetry, even if only approximate, is a highly 
desirable feature to explain the smallness of the Higgs boson one-loop correction to relation 
(p^). While the logarith mic contribution of the SM Higgs boson with a mass around the 
electroweak scale is compatible with all current precision data, any additional quadratic 
contribution is severely constrained. 

In Ref. 0, a general expression for the SU{2)l x SU{2)r chiral symmetry playing this 
important "custodial" role was presented. A phase arbitrariness in the definition of this 
symmetry has also been emphasised, in conjunction with the generic definition of the CP 
symmetry. As a result two physically different custodial limits for a CP conserving Higgs 
potential were found, characterised by the mass degeneracy of the charged Higgs pair 
with either a pseudoscalar state or a scalar state . 

The present work describes the rather unconventional phenomenology associated with 
the second implementation of the custodial symmetry, which we refer to as "twisted". 
As already mentioned in Ref. 0, this scenario gives rise to experimental signatures that 
are qualitatively very different from those usually considered in models like, for example, 
the MSSM. The main difference is the possibility for the pseudoscalar to be as light 
as a few tens of GeV with the scalar heavier than the SM expectation. This opens 
new perspectives for unusual decays like — > A^A^,H^H^,H~^H~, — > M^^^'^ and 

— > ZA^ that can become dominant and drastically affect our discovery strategies. Note 




(1.1) 
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that some of these decay modes may also appear in several BSM contexts such as a small 
tan/3 MSSM scenario (see Ref. |]), the Next to Minimal Supersymmetric SM (NMSSM, 
see Ref. [^) or the CP violating MSSM. However, current analyses of the associated signals 
appear to be limited to very specific regions of the parameter space. 

The paper is organised as follows. First, we implement a minimal version of the twisted 
scenario proposed in Ref. ||^] and discuss the possible Yukawa couplings. Second, various 
theoretical, indirect and direct constraints are examined in order to restrict the model 
parameter space and to define a pertinent set of benchmark points. Finally, production and 
decay modes are systematically considered and a few interesting experimental signatures 
identified. Three particular channels that clearly illustrate the LHC discovery potential in 
this scenario are then considered in more detail. 



2. A minimal model 



Consider a 2HDM based on the SU{2)l doublets (pi and (j)2 with the same hypercharge 
Y = +1. If one imposes only gauge invariance, the most general renormalisable scalar 
potential reads, in standard notation, 

Vi(pi, 02) = -ml4<Pi - mlct>\(t>2 - {ml4>\'p2 + h.c.) + y (</>I<Ai)^ 



+^{4>lh? + A3(0l0i)(4'/'2) + M{4>%){4>l<l^i) 



+ 



A. 



2 



(2.1) 



where m\ 2 and Ai^2,3,4 are real parameters while mg and As^g,? are a priori complex. 

This potential, though very general, has two important drawbacks. First, it contains 
a large number of free parameters that challenges any complete in-depth analysis of its 
phenomenological consequences. Second, in the absence of additional global symmetries 
like those accidentally preserved by the Standard Model scalar sector, large corrections 
to various tightly constrained precision parameters are expected. Our aim is to address 
both issues by restricting ( |2.1| ) by means of well-motivated global symmetries, yet leaving 
enough room for unconventional experimental signatures. 

Let us first assume that the electromagnetic gauge symmetry is preserved, i.e., that the 
vacuum expectation values (v.e.v.) of (j)i and 02 are aligned in the SU{2)l space in such 
a way that a single SU{2)l gauge transformation suffices to rotate them to the neutral 
components.^ After a suitable SU{2) x C/(l)y transformation, one has 



and (02) = ^ .J (2.2) 





with vi and V2 two real parameters such that vf + V2^v'^ = {■\/2Gf) ^ and V2/V1 = tan (3. 
An important feature of the 2HDM is the freedom to redefine the two scalar fields 0i and 



^This phenomenologically motivated assumption is in fact rigorously justified in the context of the 
restricted models to be considered in the following (e.g., see jl^, 
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^2 using arbitrary U{2) transformations acting in the "flavour" space, i.e., 

, U G U{2) . (2.3) 




Transformation (|2.3| ) leaves the canonicahy normahzed, gauge-covariant kinetic energy 



terms invariant. This notion of basis invariance has been emphasized in |12| and considered 



in great detail more recently in |13] and |14, 15]. From now on, let us take advantage of 



this property to select one of the Higgs basis where only one of the two Higgs fleld acquires 
a non-zero v.e.v.: 

vi = V and V2 = ■ (2-4) 

Note that the Higgs basis is not univocally defined since the reparametrization 02 e*"(/>2 
leaves the condition (^^) invariant. As illustrated in the following, this (set of) basis is 
particularly convenient since (pi contains all the SM would-be Goldstone fields, such that 
any generalization of SM global symmetries is almost straightforward. 

In the Higgs basis, the most general definition of the CP symmetry simply reads (e.g., 
see Hi): 

{CP)Mt,x){CP)^ = </>!*(*, -x) 

{CP)^2 {t, x){CP)^ = e'U2*{t, -x) (2.5) 

and displays a single arbitrary phase. For the sake of simplicity, and without any loss of 
generality, we take advantage of the (j)2 phase freedom to conventionally fix the value of 
the CP phase at one, i.e., 5 = 0. If one imposes CP invariance on the bosonic sector of 



the theory, then relations (2^) restrict all the scalar potential parameters to be real. The 
resulting potential is thus described by only ten free parameters, i.e., four less than the 
full one given in Eq. ( ^.1| ) since m| and As^gj fire now taken real. This restricted potential 
is denoted V\q^ in the following. 

In the Higgs basis, the most general custodial symmetry transformation for a two- 
Higgs-doublet model also takes the simple form 

Ml ^ UlMiUI , M2 ^ UlM2V^ (2.6) 

where Mi corresponds to the usual [1/2, 1/2] representation of (pi, i.e., 

M = ( ] (27) 

The right transformation V/j for M2 in ( |2.6D is not necessarily fixed by the right trans- 
formation Ur for Ml. Only SU{2)l x U{\)y is a local symmetry of the Lagrangian. So, 
we still have the freedom to choose Vr = X'^UrX if the two- by- two unitary matrix X 
commutes with exp(iT^), where is the diagonal generator of the global SU{2)r. We 
thus fix X = exp(i7r^) where 7 is an arbitrary angle. 

Now the crucial point is that we have already used the reparametrization invariance 
of 02 to fix the CP-phase. As an all-important consequence, the custodial angle 7 is a 
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truthfully free parameter. Depending on its value, invariance of the general CP invariant 
potential V^^^ under the custodial symmetry can be obtained for different limits of the 
potential parameters: 

1. If 7 = 0, the CP-even combinations (/'{(/'i, 4'24'2 and {(j)\(j)2 + 4>2(t)i) are invariant under 
( |2.6D . Invariance of the CP conserving potential V\l^ under this custodial symmetry 
is restored in the limit A4 = A5 such that 

+ Y(^lhf + X3{4^i)i4h) + y (402 + cpWf 
+ (0l02 + 4(Pi) [A6(0l0i) + A7(402)] . (2.8) 
This gives the usual realization of the custodial symmetry in the 2HDM. 



2. If 7 = vr, the CP-even combinations (t>\(f>i, 02 "^2 and the CP-odd combination 
0201) are invariant under ( |2.6| ). Invariance of the potential V^^^ under this "twisted" 
custodial symmetry is restored in the limit A4 = — A5 and m| = Ae = Ay = 0, giving 

F^^^*^'l(0i,02) = -ml4<P, - mUl^2 + y(0l0i)' + y(0^02)' 

+A3(0l0i)(0^02) + y (0l02 - 0^0i)'. (2.9) 

3. If 7 7^ and vr, the only combinations of 0i and 02 that are, at the same time, 
CP-eigenstates and invariant under ( p.6| ) are 0|0i and 0202. The resulting, highly 
symmetric potential reads then 

n^^"(0i,02) = -mf0t0i-mi0t</>2 + y(0l0i)' + y(0^02)' 

+A3(0l0i)(0|02). (2.10) 

In fact this potential is so constrained that the SU{2) l representation of 02 is not even 
fixed by the interaction term and any representation (singlet or higher) is allowed. 

Out of the three cases above, only the second one is novel and from now on, we focus 
exclusively on it. This "twisted" scenario, characterized by a relatively low number of free 
parameters (at least compared to the usual case), yet displays the possibility of interesting 
phenomenology due to the non-trivial mixed term (0^02 — 020i)^ in Eq. ( |2.9| ). 

The last symmetry to be considered in the twisted scenario is the Z2 symmetry 

01 ^ 01 , 02 ^ -02 , (2.11) 

known (assuming suitable Z2-charges for the SM fermions) to be a simple, yet elegant, way 
to suppress tree-level FCNC in any 2IIDM |^]. The potential (|2.9| ) is naturally invariant 
under this symmetry without further assumptions. However, if such a Z2 symmetry is only 
manifest in a Higgs basis, all fermions are forced to be Z2-even in order to couple to 0i 
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and get non vanishing mass terms. If such is the case, they cannot couple to (j)2- This 
provides a natural frame for dark matter (e.g., see [|l^]) restricting, however, considerably 
the possible phenomenological signatures at colliders. 

This apparent limitation can be circumvented if the SO (2) rotation of angle /? = 
arctan(t;2/vi) required to go from a generic basis (where both v.e.v. are real and non zero) 
to the Higgs basis (2.4) is promoted to a symmetry of the potential. In this case, a ^2 



symmetry defined in the Higgs basis would remain manifest in any related basis, and vice 
versa. Since the only matrices to commute with the generator of the SO {2) symmetry, 
i.e., the second Pauli matrix T2, are the identity matrix and T2 itself, the two invariants 
are (^j^i + 02'/'2) and (</>|</>2 — (p^'t'i)- Imposing invariance of the quartic part of potential 
( |2.9| ) under this 5*0(2) (a softly-broken Z2 symmetry being sufficient to ensure the absence 
of large FCNC effects) reduces the total number of parameters from six to four: 

Vr'" </'2) = -ml4<l>,-ml4cl>2 + ^{4h+4<l>2f 

+^i4^2-4<Pif- (2-12) 

Our forthcoming study of unconventional phenomenology is based on this minimal potential 
where, for convenience, we have introduced the symmetric (Ag = Ai = A2 = A3) and 
antisymmetric {Xas = A5 = — A4) notation for the coefficients of the corresponding quartic 
terms. 

It is straightforward to determine the physical spectrum associated with the po- 
tential ( |2.12| ) in the Higgs basis. First, it contains a CP-even SM-like Higgs boson 
= V2(Re{4)i) — v/\/2) with squared mass 

mlo = Xsv'^ = 2m\ . (2.13) 

Second, it displays a pair of charged Higgs bosons and a CP-even scalar F° = -y2Re(02), 
forming a triplet under the custodial symmetry with 

2—2 2 22 /0 1/l^ 

m-j. = vn^^ = rrij^o = w-i — ■ (2.14) 

Last but not least, it allows the pseudoscalar state = -v/2Im(i?!)2) to be singlet under the 
twisted custodial symmetry with 

m\o = rn\ - Xasv^ ■ (2-15) 

Since hP remains the only massive Higgs boson in the limit of an exact 50(8) symmetry 
acting on the real doublet components [Xas 0, mi — > 7712), one may thus expect the 
following inverted hierarchy: 

TJiy^o < mfjO [j± < niyfi if Xas > 0. (2-16) 



As seen from Eq. ( 2.14 ), the main consequence of the twisted custodial symmetry 
is a mass degeneracy between the pair of charged Higgs bosons and the scalar state 
H^. This contrasts with the usual approach that links the custodial symmetry to a mass 
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degeneracy between the pair of charged Higgs bosons and the pseudoscalar state . 
The existence of a genuine twisted scenario in the 2HDM, and more generically the in- 
terpretation of the possible interplay between the custodial and CP symmetries in terms 
of aligned/misaligned phases, Eqs. (2.8)/(2.9), is a non trivial observation. Previous phe- 



nomenological studies of the p ~ 1 constraint had already noted a vanishing contribution in 
the 'yn'\j± ~ "T-^o limit or more recently and [l^) but none of them did interpret 
this property in terms of symmetries. The interplay between the CP and custodial sym- 



metry in 2HDM has also been considered in [20|, but the second CP conserving custodial 



scenario discussed there is in fact physically indistinguishable from the "usual" case since 
can be redefined as CP odd. Here, the usual and twisted cases clearly lead to differ- 
ent physical situations, as will be shown explicitly by considering, for example, quantum 
corrections to the T parameter (see Fig. ^) or to the muon anomalous moment (Sec. 3.2.5| ). 



Regarding Yukawa couplings, the presence of a softly-broken Z2 symmetry allows to 
define type I and type II models |^], with a very different phenomenology. Type I models 
correspond to the case where only one Higgs doublet couples to all fermions in the generic 
basis. This can be enforced by assuming that all left and right handed fermions are even 
under the softly-broken Z2 symmetry (2.11). The Yukawa interaction terms in the Higgs 
basis then read 

/2 ^ 

£f = QLMdicpi - tan/3 (/-ajd^j QlM„(0i - tan(3(P2)uR + h.c. , (2.17) 

V V 

with 01 = iT24'\. In terms of physical states, the couplings of to fermions are identical to 
those observed in the SM, i.e., directly proportional to the fermion mass matrix M, while 
those of H^, and are simply rescaled by a factor^ tan/3. If tan/3 — > 0, (/!>2 decouples 
from the fermionic sector and provides a viable dark matter candidate ||l6|| . 

In type II models, one of the two Higgs doublet couples to the down type quarks and 
to the charged leptons, while the other one only couples to the up type quarks in the 
generic basis. This is easily achieved by assuming all fermionic fields to be even under the 
softly-broken Z2 symmetry ( p. Ill ) except for the right handed up type quarks fields which 
are odd. The Yukawa interaction terms in the Higgs basis now read 

/2 ^ 

Cif = -—QLMd{(t)i - tan/3 (^2)dR - —QlMu{^i + cot /? (^2)ur + h.c. . (2.18) 

V V 

In terms of physical states, the couplings of to fermions are still identical to those 
observed in the SM, but those of H'^, A^ and are now rescaled by a factor tan /3 or 
cot /3 depending upon whether the fermionic current involves down or up type fermions. 

The minimal two-Higgs-doublet model (M2HDM) considered in this work is defined 
by Eqs. (|2.12| ), ( 2.17 ) and ( |2.18|) . It has only four free parameters: mj-fi, uit and m^o for 



the scalar potential, and tan fi for the Yukawa interactions. Inspired at first by a twisted 
custodial symmetry, one can also find a generic 2HDM whose parameters in the potential 
are adjusted in such a way to implement the inverted mass spectrum ( p. 16 ). 



^Depending on which doublet is conventionally chosen to couple to all fermions in the generic basis, this 
scale factor can also be cot /3 as often considered in the literature. 
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3. Parameter constraints 



In order to reduce the parameter space, let us now review the theoretical, indirect and 
direct constraints at our disposal. 

3.1 Theoretical constraints 

The first type of relevant theoretical constraints are the vacuum stability conditions. They 
come from the requirement of a positive potential for large classical values of the fields in 
any direction in the (j)2) plane. These constraints can be obtained by considering only 
the quartic terms of the potential p^ , In the context of the minimal potential defined 
in Eq. (|2.12D , they read 

A5>0 and A5 > Aa5 , (3.1) 
or equivalently, in terms of the physical masses, 

m^o > and m^o > tti^ — m\o . (3.2) 

In fact, the squared mass of all the (pseudo)scalars are, by definition, positive if the v.e.v. 
correspond to a minimum of the potential ( p. 121) . 

The second set of relevant constraints comes from the requirement of unitarity for all 
the possible scattering processes involving the new scalar particles. They have been worked 
out for both CP conserving and CP violating potentials ||2^, |2^ and can be advantageously 
summarized as 

|A^2^3|<87r, (3.3) 

where | Ay^^g | are the eigenvalues of the high energy scattering matrix for different quantum 
numbers of the initial state: total hypercharge Y, weak isospin and Z2 parity. In the 
M2IIDM (see Eq. ( |2.12| )), the relevant contributions are 

yYeven = -|- 

Air = 3A5 ± (2As - Xas) 

Aoo'' = As - 2Xas ± 3Aas • (3.4) 

The two sign possibilities appearing in these expressions correspond to the scattering of 
different initial states with the same quantum numbers. Using relations ( 2.13| ) and ( 2.15| ), 



constraints ( p.3| ) restrict the possible values of the scalar masses. In particular, in the limit 
where all scalar masses are small except for one, say ms, one has the upper bound 

m5< 550GeV. (3.5) 

When the Higgs masses are non negligible, the unitarity requirement may help to restrict, 
for example, the allowed region in the (771^0,771^) plane for given values of 777/jO, as displayed 
in Fig. |^(a). A naive estimate of the region compatible with the perturbative approach 
is shown in Fig. |^(b). In this region, all three- and four- scalar vertices are bounded by 
47r, such that the effective parameter of perturbation theory is smaller than one. Let us 
emphasize that more sophisticated treatments of the perturbativity constraint applied on 
similar models (SM and MSSM) only lead to minor corrections with respect to this naive 
approximation (e.g., see p5| ). 
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(a) Unitarity (b) Perturbativity 

Figure 1: Unitarity and perturbativity constraints in tlie (m^o, m^) plane for the M2HDM. Dotted 
red lines are limits for my^o = 120 GeV, dashed green lines for m/jO = 300 GeV and plain blue lines 
for m/jO = 500 GeV. The allowed regions lie between these lines. 



3.2 Precision measurement constraints 
3.2.1 Electroweak precision parameters 



The total contribution of the new scalar states to the T parameter [26| in the context of the 



multi-Higgs-doublet model has been computed recently in Ref. [27|. Here, we focus on the 
correction AT from the CP conserving potential V^J^. Those are well-known in the limit 
where all scalar squared masses are larger than |17, 28 1. If one of the scalar is lighter 
than mz, however, the exact expression obtained in Ref. [^] and reported in Ref. pO| ] is 
more suited. 

The results for AT as a function of m^± are shown in Fig. With different values 
of (/5 — a), where a is the mixing angle between the two CP-even eigenstates in and 
fixed values for all other scalar masses. Fig. ^ displays two cases where AT is close to zero: 

1. The solution mj|^± m^o, independently of the value of (/? — a) and the masses of 
the other scalars; 

2. A continuum of solutions ranging from mu± ~ mfjo when (5— a = vr/2, to mu± ~ rriyfl 
when /? — a = 0. 

The first possibility corresponds to the usual custodial scenario leading to the potential 
( |2.8| ), while the second one corresponds to the "twisted" case with the potential ( |2.9| ). In 
fact, the two solutions mu± = ra^fi and mu± = mjjo in the second case are related through 
a 7r/2 shift in a, i.e., a renaming of and H^. The continuum of solutions between these 
two mass degeneracies corresponds to the situation where the neutral state belonging to 
the custodial triplet is not a mass eigenstate, but a mixture of and H^. 
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AT 




Figure 2: The AT correction with respect to the charged Higgs pair mass (in GeV) from the CP 
conserving potential in a 2HDM, for (/? — a) = (dotted red), 7r/4 (dashed green) and tt/2 (sohd 
blue). The other scalar masses are fixed to nifio — 150 GeV, m/^n = 400 GeV and m^o = 1 TeV. 
The thin horizontal black lines shows the 2a experimental limits on AT (assuming AS ~ 0) 



An interesting possibility in the framework of the M2HDM based on potential ( 2.12| ) 



arises when the pseudoscalar state is light while all the other scalars are heavier (> 
100 GeV). In this case, a moderate deviation from the degeneracy 'mfj± = rrifjo could 
compensate the large logarithmic contributions involving m/jO. This can be seen directly 
from the first order analytical approximation for AT in the mH± ~ nifjo region ||lJ 

2 2 



- 2 2^ X < COt^ Ow 

Ibvrm^ cos^ Uw 



-3n.^(log^ + -i-log^ + l)^ . (3.6) 



The amount of breaking required for a vanishing T as a function of the Higgs mass is 
shown in Fig. ^ Using the effective potential approach [^], we estimated the one-loop 
corrections to the tree- level relations ( 2.14| ) in the context of the M2HDM. Non- vanishing 



contributions can come only from terms that explicitly break the custodial symmetry, i.e., 
the gauge and Yukawa interactions. The main contribution arising from the top (and 
bottom) quark loops leads to a relative mass difference {mjj+ — rrijjo ) /m jjo smaller than 
1% for both type I and type II models. The remaining mass difference required to satisfy 
the T parameter constraint (see Fig. ^) should thus either be explained by the presence of a 
small custodial breaking term in the scalar potential, like for example e{(l)\4>2+4>24'i)'^j or by 
the presence of other BSM contributions {e.g., heavy fields leading to additional radiative 
corrections). Yet, those effects would have a negligible impact on the phenomenology 
beyond AT. 

Contrary to the T parameter, the S and U parameters only depend logarithmically 
on the mass of the new scalar particles. The exact one-loop contribution from V^J^ has 



been evaluated in Refs. [29, The numerical results for these contributions show that 



scenarios with a light pseudoscalar and a heavy degenerate triplet are favoured by both 
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Figure 3: Soft relative corrections to mj/± = mno in the M2HDM that give AT = as a function 
of the SM-hke Higgs boson mass m^o (in GeV). The dotted red, dashed green and sohd blue hnes 
correspond to mjLfo=200, 300 and 400 GeV, respectively. The mass is fixed at fOO GeV but does 
not affect the results if m^o ^ »ti^o, m^. 



S and U parameters. In this region of the parameter space, the 2HDM contribution 
has an opposite sign compare to the SM one. For sufficiently small pseudoscalar masses 
(and large triplet masses), this extra contribution could even partially compensate a large 
positive contribution to the S parameter induced by a heavy (« 300 GeV) SM Higgs. A 
contrario, scenarios with a very heavy pseudoscalar and a light triplet are disfavoured or 
even excluded (depending on the actual value of the T parameter) by the experimental 
upper bounds on S. 

3.2.2 Rare B mesons decays 

The b ^ branching rate measurement [^] is known for its stringent bounds on charged 
Higgs boson masses. The leading contribution of a charged Higgs boson to this loop- 



induced decay has been derived in Refs. |34, 35]. As expected, the small tan/3 region 



for the type I 2HDM ( p. 17 ) is left unconstrained since, in this case, the charged Higgs 
bosons decouple from the fermions. For larger values of tan f3, only a very narrow region 
of the parameter space survives the constraint. In type II 2HDM ( p. 18 ), the leading order 



prediction for the lower bound on the charged Higgs mass (> 500 GeV at 95% CL) is 
essentially independent of tan /?, for tan /? > 2, due to the cancellation between the tan /? 
and the cot P contributions. 

However, this leading order prediction of the b ^ sj branching ratio suffers from large 
uncertainties which can be partially reduced by taking into account higher order QCD 



corrections. For the 2HDM, an estimation of the NLO corrections in |3(:] shows a sizeable 



effect which is sufficient to drastically reduce the lower bound obtained at LO for the 



mass. The current status for a type II model is summarised in Ref. |37]. The 95% (99%) 
lower bound on mjj± is around 295 (230) GeV and stays practically constant down to 
tan/3 ~ 2. Experimental results may even be interpreted as favouring a charged Higgs 
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mass around 650 GeV. As already mentioned, type I scenarios at low tan (3 values are 
close to the decoupling limit, such that the new physics corrections are generally small in 
magnitude (but of opposite sign compared to type II) and the NLO effects are not relevant 
(e.g., see Ref. for a discussion). At higher tan /? values, the strong coupling regime 
is quickly reached and even the NLO prediction may not be well-behaved (i.e., very scale 
dependent) in this region. In the following, we only consider the tan /3 < 0.5 region and 
discard the possibility of a high tan /? value. 

If the charged Higgs boson couples strongly enough to the be or bu quark currents, and 
at the same time to the rlv leptonic current (i.e., in type II scenarios with large tan /?), the 
B DrVr and B rvr branching ratios could be affected. The normalised branching 
ratio 

measured by the BaBar collaboration [^] sets a 95% C.L. bound on the type II 2HDM 
parameter space, namely 

tan (3 < 0.23 GeV"^ x m^^± . (3.8) 

Taking into account the lower bound on mjj± from the b ^ process, this last constraint 
can only be relevant in the very high tan [3 region (tan (3 > 70) which should anyway be 
discarded by the requirement of perturbativity for the Yukawa couplings. 

A similar bound can be obtained from the rare process B — > tVt. The experimental 
result from Belle [^] can be compared to the best SM prediction as follows 

BR(B TVr) , . 

^ ' -1.13 ±0.44 (3.9) 



BR(5 ^ TVr)sM 

where only the experimental error is considered. A normalization with respect to the 
2HDM type II model (for a review with Minimal Flavour Violation, see [^]) gives 

tan /?< 0.13 GeV"^ X mj^± . (3.10) 

Assuming m/^± > 300, this constraint may significantly restrict the tan/3 > 40 region. 
Let us however emphasise that this bound should be considered with some caution: the 



B TVr signal "evidence" in the BaBar experiment |43| is still statistically lower than 
that obtained by the Belle collaboration, and the theoretical uncertainties associated with 
the lattice estimate of could be underestimated. 

3.2.3 The Bq-^ mixing 

The virtual effects of the charged Higgs bosons on the — Bd oscillations have been 



described at leading order in Ref. ||4J. The resulting bounds on the charged Higgs mass, 
with respect to tan P, are visible on Fig. |^. Similarly to 6 — > 57, and in order to avoid 
discussing the actual choice of the SM parameters values entering the Bq — Bq mixing, we 



normalise it to recover the most recent SM prediction from lattice simulations [45| 

Am|^ = 0.69 ± 0.15 ps~^ (3.11) 
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tan jS tan (3 

(a) Type I (b) Type II 

Figure 4: Bq — Bq mixing bounds on the charged Higgs mass (in GeV) with respect to tan/? in 
the 2HDM with (a) type I and (b) type II scenarios, at two (dashed green) and three (soHd blue) 
standard deviations. There are no la limits on these figures due to the slight discrepancy between 
the SM prediction and the current experimental measurement. Predictions are at leading order. 



to be compared to the experimental world average [46| 



Am'^^P = 0.507 ± 0.005 ps"^ . (3.12) 

Consequently, the error on the resulting constraint is clearly dominated by the theoretical 
uncertainty in ( 3.11] ). For type I scenarios, the tan/3 < 0.5 — 1 constraint obtained on the 



whole range for the charged Higgs mass is similar to the one obtained from the b ^ 
process. In the case of type II 2HDM, the tan P < 2 region is excluded at more than 95% 
CL, almost independently of mfj±. Like in the 6 ^ 57 case, the inclusion of 0{as) QCD 
corrections reduces the sensitivity of Ams to charged Higgs contributions and slightly 
weakens the above constraints 

3.2.4 The Zbb vertex 

Loop corrections involving new charged and neutral scalars may give sizable contributions 
to Rb, the observable hadronic branching ratio of Z bosons to 66, and to A/,, the 6-quark 



asymmetry. These corrections have been derived in great detail in Ref. |48] for the CP 
conserving 2HDM. Corrections to Rb and Ab are expected to have similar magnitudes. 
However, the high experimental precision associated with the Rb measurement makes it 
much more discriminating than Ab on the whole parameter space, as shown by an explicit 
numerical analysis. 

In type I models, the only relevant contribution is the charged Higgs correction to gi 
since all the other ones are suppressed by nib/ v. However, the typical bound extracted 
in this case, namely tan/3 < 1, is not competitive with the one coming from the Bq — Bq 
measurement. The situation is more interesting in type II models where the neutral scalar 
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contributions are potentially sizable. In the type II M2IIDM, this restricts the allowed 
region in the {mj^o,mfjo) plane, as illustrated on Fig. |5|. Since the {H^,H^) triplet is 




100 200 300 400 500 100 200 300 400 500 

(a)tan/3 = 30 (b) tan/3 = 50 



Figure 5: Constraints at la (dotted red), 2a (dashed green) and 3a (solid blue) on the type II 
M2HDM parameter space from the Rb measurement, for two different values of tan/3. 



forced to be rather heavy by the B physics constraints previously reviewed, the mass 
of the pseudoscalar is bounded from below. For tan /3 = 50 and m/^± > 300 GeV, 
for example, the bound is approximatively m^o > 60 GeV at 95% CL. This bound is 
somewhat lower than the one obtained in Ref. [^8| (m^o > 100 GeV) in the days when the 
experimental value of Rb was sizably larger, leaving less room for negative contributions 
from new physics. 

3.2.5 The muon anomalous magnetic moment 

The anomalous magnetic moment of the muon is known to be a rather interesting probe for 
BSM physics. The SM theoretical prediction, incorporating the e"'"e~ — > vrvr data obtained 
by CMD-2, KLOE and SND, gives |4|] 

^th,SM ^ (11659180.4 ± 5.1) x 10"^° (3.13) 

to be compared to the experimental measurement from the Brookhaven experiment |5^] 

a^^P = (11659208.0 ± 6.3) x 10"^° . (3.14) 

The 3.4(7 deviation between these two values may be optimistically interpreted as a signal 
of new physics, or, at least, as a valuable constraint. Due to the smallness of the muon 
mass, this constraint is only relevant for Higgs physics when the coupling to the leptons is 
increased compared to its SM value. In the M2IIDM considered here, this corresponds to 
the type II scenario for Yukawa couplings, for which the one- and two-loop contributions 
are reported in Ref. |51|. 
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At one-loop, the scalar contribution ajj is positive whereas the pseudoscalar and the 
charged Higgs boson give negative contributions. Each contribution reaches its extremum 
at small masses and vanishes like m^/m| log(m|/m^) at large masses. The absolute mag- 
nitude of each type of contribution is shown on Fig. |6|(a) for tan /? = 1. The total one-loop 




(a) (b) 

Figure 6: (a) Absolute value of the one- loop contribution to from a neutral scalar (dotted 
red), a pseudoscalar A" (dashed green) and a charged Higgs boson (solid blue), if tan/3 = 1. The 
contribution is positive while the A" , are negative, (b) Same for the two-loop contributions 
from a neutral scalar (dotted red) and a pseudoscalar A" (dashed green). At two- loop, the 
neutral scalar contribution is negative while the pseudoscalar one is positive. Only the 6, r and fi 
fermion loops are included. 

correction is dominated by the contributions of neutral Higgses for masses above 0.2 GeV. 
Solving at the one-loop level the theory/experiment discrepancy in the M2HDM with 
a moderate tan/3 would require a very light (< 10 GeV) scalar . The situation changes 
when considering the dominant two-loop contributions. The contribution is now neg- 
ative whereas the one is positive. The total two-loop contributions can be seen in 
Fig. ^(b). By comparing with Fig. ^(a), one finds that these two- loops corrections are 
dominant for ms ^ 10 GeV (they cancel against the one-loop part for ms ~ 5 GeV). 
Due to the opposite sign of the scalar and pseudoscalar contributions, solving the the- 
ory/experiment discrepancy in the M2HDM with a moderate tan/? in this region would 
require a light pseudoscalar (20 < m^o < 100 GeV). If tan/3 ~ 30, a very good agreement 
can even be reached for m^o ~ 20 GeV. A larger tan j3 value would require a heavier 
pseudoscalar and vice versa. In all cases, corrections to clearly differ when considering 
either scalar or pseudoscalar particles. The a^ measurement can thus potentially help to 
physically disentangle the usual and twisted custodial scenarios introduced in Section |2[ 

3.3 Collider constraints 
3.3.1 LEP 

LEP searches for a SM Higgs boson in the standard decay modes hP bb and — > r+r", 
have provided an exclusion at a 95% confidence level for a Standard Model Higgs boson 
with a mass lower than 114.4 GeV ||5^. In the context of models with an extended scalar 
sector, may also decay, possibly dominantly, in a pair of lighter Higgs bosons. Since 
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such a light object has escaped the LEP searches, it must either have a reduced couphng 
to ZZ or unusual decay properties. In the context of the M2HDM, both and could 
satisfy this requirement thanks to their vanishing couplings to Z boson pairs. However, 
here is nearly degenerate in mass with H^. Since the possibility of a light (< 100 GeV) 
charged Higgs boson is strongly disfavoured by both direct and indirect measurements, we 
then focus on the light hypothesis in the following. 

In the forthcoming phenomenological analyses, we conservatively consider the usual SM 
bound wiyfi > 114.4 GeV. However, it should be mentioned that this bound could be slightly 
lowered if the — > bb branching ratio was reduced due to the presence of the A^A^ 
decay mode. Constraints on this decay mode from LEP data have been considered, for 
example, in the framework of the Next-to-Minimal Supersymmetric Standard Model (see 
|[5^ and reference therein). If my^o > 2mh, it has been first thought that this decay could 



explain the simultaneous excesses observed in the Z2b and ZAb |54] final states by 
adjusting m/,o, BR(/i° BR(/i° bb) and BR(^° bb). However, the ZAb excess 

tends to favour slightly higher masses for (in the 105-110 GeV region) compared to the 
main Z2b excess (in the 100 GeV region), decreasing the significance of a global fit. 

The model independent searches for the e~^e~ — > Z* ^ H^A^ pair production process 
at LEP (in the 46, 262r and 4t channels) put the tightest constraint on the M2HDM 
mass spectrum, and in particular on the type I model where indirect constraints from B 
physics still allow for a relatively light triplet {H^,H^). Taking into account the values 
of BR(i?o, AO ^ bb) and BR{H^,A^ r+r"), the final result from [||] can be roughly 
summarized as the bound m^o + mjjo > 170 GeV for m^o mjjo. If AO is very light, 
namely m^o < 30 GeV, a slight loss of efficiency in the 46 and 262r analysis allows for 
rriHO > 130 GeV. 

If niHO > mz + rriAo, the scalar could also decay into ZA^, thus reducing dra- 
matically the branching ratio of T^T~,bb. The final state signatures associated 
with this possibility (namely Z4b, Z262r and Z4t, see also Ref. for a discussion in the 
small tan (3 MSSM) would mimic the A^A^ process described previously, but with a 
different kinematical structure. In this case, the decay products of the Z boson and the two 
softer 6-jets should be used to reconstruct the mass. In the absence of any dedicated 
experimental study for this open possibility, we adopt a conservative approach and restrict 
the model using the (mj:^o,m^o) limits already mentioned. 

Another strong constraint in the low m^o region for type II models can be deduced 
from searches for the Yukawa-induced process e+e" — > 66ylO with A^ r^r^^bb. The 
result for each mode is available in Ref. |55|. Due to the reduced branching ratio for 
AP t'^t~, the A^ — > 66 mode is the most restrictive one in the rrij^o > 2mb region. The 
limit is tan /? < 20 for m^o ~ 10 GeV and quickly becomes less relevant for higher masses 
due to the smaller production cross section. The loop decay Z A^'y (through a quark 
or a lepton loop) could in principle also be used to constrain the A^ mass. However, an 
extensive analysis of this channel [^6| has shown that the LEP measurement sensitivity (of 
order 10~^ for the associated branching ratio) was not sufficient to put a lower bound on 
m^o tighter that the one obtained from the Yukawa-induced process. 

Finally, results for charged Higgs boson searches at LEP in the general 2HDM (type 
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I and type II) by LEP collaborations |57, motivated by the possibility for three-body 
decays of Higgs bosons reviewed in |59], are also considered. In addition to the usual 
fermionic decays — > t'^i'j- and H'^ cs, the possibility for charged Higgs boson pair 
production (e+e~ — > Z* ^ H^H^) with the decay into has also been taken into 

account. The existence of a charged Higgs boson with mass lower than 76.7 GeV (type 
I) or 74.4 GeV (type II) is excluded at the 95% CL, for a wide domain of the parameter 
space. 



3.3.2 Tevatron 

To date, the total integrated luminosity collected at the Tevatron collider is not sufficient 
yet to exclude a Standard Model Higgs boson at masses higher than the LEPII limit, except 
for a narrow window around 165 GeV [^]. We thus focus here on the — > j4'^A'^ exotic 
decay, with — > bb,T~^T~. 

The case of direct production of through its effective coupling to gluons, followed 



by the decay — > A^A^ — > 46 has been covered in Ref. |61|. Unfortunately, the Ab 
QCD background overwhelms the signal and a discovery could only be achieved if the 
production would be enhanced by at least one order of magnitude. The same process where 
one of the A^ decays into t^t^ instead of bb could provide a cleaner signature, at the price 
of reducing the rate due to the small A^ t~^t~ branching ratio. 

Since the associated production of hP with a vector boson V = W, Z \s the second 
largest production mechanism at the Tevatron, it is also natural to consider the exotic 



decay hP — > A^A^ in this context. Detailed studies (see Refs. [62], [33|, |64| and the 
references therein for an overview) of both the V2b2T and VAb final states have proven the 
potential for such signature. However, the associated statistical significance is too small to 
constrain the model parameters with the currently available integrated luminosity. 

Searches for the and A^ bosons at the Tevatron experiment take place for pro- 
duction in association with b quarks, or in gluon fusion, and decays to bb and t~^t~ final 
states. Since most analyses are oriented towards the MSSM Higgs bosons discovery, they 
focus mainly on the ms > 70 GeV mass region. Exclusion regions for Run I |^5[ and Run 
II [^] can be (very) conservatively translated as a constant tan /? < 35 bound on the whole 
^Ao,HO > 70 GeV mass range. 

Another interesting possibility to produce a light A^ boson at the Tevatron (and pos- 
sibly at the LHC) is the charged Higgs associated production pp W^* — > H^A^, where 
the charged Higgs may further decay into Vl^^A'', which has been proposed in Ref. |67] and 



discussed more recently in Ref. [|68| in the context of the NMSSM. In the most favourable 
mass scenarios (m^o < 20 GeV, mfj± < 90 GeV), the associated cross sections can be 
larger than 500 fb and a signal could be detected in the Wib final state. 

Finally, if mjj± < rn-t, a top quark could have a sizable branching ratio to H~^b. This 
possibility has been considered by the Tevatron experiments for various charged Higgs decay 
hypothesis. The result is a 95% CL upper bound on the t H^b branching ratio, i.e.. 



0.2 for the whole mass range [69| 
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Bq — Bq mixing 


tan /? < 35 


B TVr and LEP hhA 46 



Table 1: Summary of the relevant constraints for the M2HDM defined by Eqs. ( |2.12| ), ( ^.17D and 
(H). 



3.4 Constraint summary and mass spectrum 

All the relevant constraints introduced in the previous sections are summarised in Tabjl|. 
Only the most stringent bounds are presented, and some of them could be strongly cor- 
related with others. These constraints have been derived with the implicit assumption 
that all BSM contributions beyond the M2HDM are negligible. They thus appear more 
as guidelines to restrict the parameter space for the phenomenological study of collider 
signatures presented in Chapter 4, than as strict bounds. 

One can foresee qualitative differences between the "twisted" scenario and other 2HDM 
realisations. The (constrained) CP-conserving MSSM Higgs phenomenology in the low 
mass region is very similar to what can be expected in the SM. So, an identification of its 
enlarged scalar sector at the LHC will have to rely on the direct detection of the heavy 
states H^, and AP , which may require a high luminosity (~ 100 fb~^, depending on 
tan (3). The same conclusion also holds for most of the "usual" scenarios considered in the 
literature, where the nearly degenerate custodial triplet {H^,A^) is forced to be relatively 
heavy by severe B meson physics constraints. 

The situation could be completely different in the M2IIDM considered here, due to 
its distinctive inverted spectrum (see Eq. ( 2.16| )). In type I scenarios, the reduced Yukawa 
coupling of the charged Higgs allows for a moderately light custodial triplet {H^,H^). 
A small mass splitting inside this triplet (and, to a lesser extent, the presence of a light 
pseudoscalar A^) is required to help a rather heavy SM Higgs (say, m^o « 300 GeV) in 
passing the electroweak precision tests. If mf^o > 2mT, many exotic scalar decays can 
be kinematically allowed {h? — > H~^H~ , hP H^H^, . . . ) giving rise to unconventional 
collider signatures. In type II scenarios, the custodial triplet mass is constrained to be 
larger due to the 6 ^ 57 bound on the charged Higgs mass. A relatively small mass 
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splitting inside this heavy triplet (of order 5%) is enough to allow for a very heavy SM 
Biggs'^ (m^o ~ 400 GeV) . The electroweak oblique parameter S together with the available 
data for the muon anomalous magnetic moment slightly favour the presence of a light 
pseudoscalar in this case but its mass is bounded from below (m^o > 30 GeV) by the LEP 
direct searches and the measurement. Even though the m^jO > 2m condition is hardly 
satisfied in this context due to unitarity and perturbativity constraints on mf^o, unusual 
decays for and remain an interesting open possibility. 

4. Discovery potential at the LHC 

This last section is dedicated to the discovery potential of type I and II M2HDMs at the 
Large Hadron Collider (LHC). First, we review the main Higgs production and decay modes 
for both type I and type II scenarios. Next, the constraints derived in Section ^ are used 
to identify two benchmark points in the allowed regions of the model parameter space, one 
for each type. Finally, after a short description of simulation and reconstruction tools, the 
prospects for three particularly relevant signal processes at the LHC are reviewed in more 
detail. 

4.1 Production and decay modes of Higgs bosons and benchmark points 
4.1.1 Type I scenarios 

The type I M2HDM is characterised by a tan /3 parameter smaller than unity. In this case 
the only relevant direct production mechanisms for H^, and involves couplings to 
top quarks and/or gauge interactions (e.g., qq ^ Z* ^ H^A^), while all the usual SM 
production modes are still relevant for the h^. 

In general, the cross sections for the production of pairs of scalar particles barely reach 
1 pb. For H^H^ and H^H^, it is even smaller due to the constraint on their mass to be 
heavier than 130 GeV. Around this mass bound, the dominant production mechanism for 
the is via the top decay t H^h due to the very large tt production cross section at 
the LHC, 

limij± > 150 GeV, the associated production of the charged Higgs boson with a single 
top quark becomes dominant compared to the t H^b decay which goes down quickly 
as mjj± gets closer to the kinematical bound mt — nib. Finally, the production of neutral 
Higgs bosons associated with the ti pair is also possible but cross sections are smaller than 
for the SM Higgs boson. 

The main decay modes of the h^, shown in Fig. |7| for two different values of the A^ 
mass, can be qualitatively very different from those observed in the SM or in the MSSM. 
In the low mass region, i.e., if m/jO < 2m^o, the primarily decays into b, r and c pairs 
(in order of importance) like in the SM. If 2m^o < m^jO < 2mw, the dominant decay 
mode is — > ^''^'^ due to the large trilinear scalar coupling. Note however that this 



'^Evcn though it appears as an interesting possibility required in the context of the M2HDM to satisfy 
Eq. (3.2), let us emphasize the presence of a heavy SM Higgs boson plays no role in the phenomenological 
signatures of type II scenarios considered in Sections 4.3 and i.4. 




(a) m^o = 30 GeV (b) m^o = 70 GeV 



Figure 7: Tree- level branching ratio of /i° into fermions (dotted lines), vector bosons (dashed lines) 
and scalars (solid lines), with respect to its mass and for two different values of m^o. 

coupling, and hence the — > branching ratio, vanishes in the narrow mass window 

around ruf^o = yj2{mj. — m^o)- For 2mw < mh^ < 2mT, the decay into a pair of gauge 
bosons dominates with BRs similar to the SM ones. Above the 2mT threshold, both the 
H^H^ and the h? — > H^H~ decays are kinematically allowed and account for about 
half of the total width. This total decay width is sizably larger than in the SM at low 
masses (around 1 GeV for m^o = 100 GeV), due to the unusual —>■ A^A^ decay mode, 
but it remains of the same order of magnitude at higher masses. 

For the other Higgs bosons, A^, and H^, the only possible decays are into SM 
particles (scaled by tan/3) if rn}j±/j^o < ^T^vF±/z + rnj^o. If the (H^) mass is above the 
ZA^ (W^A^) mass threshold, the branching ratio (BR) into these particles become close 
to unity. For and H^, the total decay widths remain at most of the order of the SM 
Higgs one at the same mass. 

4.1.2 Type II scenarios 

All production mechanisms for h^, and those involving only gauge or scalar interactions 
for H^, A^ and H^, are identical to those described above for type I scenarios. The main 
difference is the enhanced coupling of the extra Higgs bosons to down-type quarks and to 
charged leptons. The direct production through b quark annihilation thus becomes by far 
the main production mechanism at the LHC for the neutral Higgs bosons and A^. 

The couplings remain the same in type II and type I models, hence the branching 
ratio patterns shown on Fig. remain valid. The and decays are still dominated 
by the ZA^ and W^A^ at high masses. However, contrary to what happens for type I 
models, the decays into bb and tb respectively are not negligible and become dominant for 
masses below 250 GeV. 

4.1.3 Benchmark points 

The identification of the unconstrained region for the parameter space of type I and type 
II scenarios drives the choice of benchmark points in the Higgs inverted mass spectrum 
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(see Eq. ( ^.161) ). The possible values for mx and uia, taking into account direct and 
indirect constraints, are summarized in Fig. ^ In this work we restrict ourselves to two 




(a) Type I (b) Type II 



Figure 8: Summary of all relevant theoretical (green), indirect (red) and direct (blue) constraints 
on the type I and type II M2HDM, in the plane [myi,my]. Gradient bounds indicate milder 
indirect constraints, like constraints associated with the S and parameters, which should not be 
considered too strictly. The purple stars indicate the two choices of benchmark points as detailed 
in the text. 



representative cases, one for each type: "BPl" for the type I and "BP2" for the type II. 
For each case, the complete set of relevant parameter values is given in the Tab. 

4.2 Event simulation, reconstruction and selection cut definitions 

Signal and background events have been simulated using the generic 2HDM model of the 



tree- level matrix-element based event generator MadGraph/MadEvent v4.4 |70, 71 1. The 



parameters of the model have been calculated using the TwoHiggsCalc calculator jTOj. 



The PDF set used is CTEQ6L [72| and the factorisation (/^f) and renormalisation (^j?) 



scales are evaluated on an event-by-event basis using relation 

^,1 = ^1 = {MI,, + Y,Pt) (4.1) 

i 

where Mmax is the larger mass among the final state particles and j runs over the visible 
particles. 

The showering/hadronisation phase, as well as the decay of unstable SM particles, are 
simulated using Pythia 6.4 [ |73[ . 

In order to take into account the efficiency of event selection under realistic experi- 
mental conditions, the fast detector simulator Delphes [74| has been used. Characteristics 
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Table 2: Parameters values for the M2HDM that define the BPl and BP2 benchmark points. The 
branching ratios relevant for the analyses presented in the following are also given. For simplicity, 
the small triplet mass splitting has been neglected in the BP2 case. 



of the simulated detector, i.e., its geometry, granularity and resolution, are close to those 
associated with the ATLAS and CMS detectors. The tracker is assumed to reconstruct 
tracks within \ri\ < 2.5 with a 100% efficiency and the calorimeters cover a pseudo-rapidity 
region up to |r/| < 3 with an electromagnetic and hadronic tower segmentation of Ar] ~ 0.1 
and A(p ~ 0.1. The energy of each quasi stable particle is summed up in the corresponding 
calorimeter tower. The resulting energy is then smeared according to resolution functions 
assigned to the electromagnetic calorimeter (EC) and the hadronic calorimeter (HC) pa- 
rameterized by: 

"^C „.o„5+2# + 2:2£ (4.2) 



E E \/E 

^^.0.05.1|. ,4.3) 



where the energy E is expressed in GeV. 

The acceptance criteria are summarized in the Tab. ^. For the lepton, we demand a 
tight isolation criterion: the number of additional tracks with Pt > I GeV (denoted ]^tracks) 
present in a cone AR = \J Arj^ + Acj)^ = 0.3 centered on the lepton track must be either 
or 1. In so doing, we accommodate cases where collinear hard leptons are produced. The 
jets are reconstructed using only the calorimeter towers, through the SISCone algorithm, as 
defined in the FastJet package |^ and implemented in Delphes. Unless stated explicitly, a 
cone size radius of 0.7 is applied for the jet algorithm. The 6-tagging efficiency is assumed 
to be 40% for all 6-jets, independently of their transverse momentum, with a fake rate 
of 1% (10%) for light (charm) jets. Finally, the total missing transverse energy is 
reconstructed using information from the calorimetric towers and muon candidates only. 

In order to avoid repetitions in the forthcoming analyses, we define here a set of cuts: 
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Final state 


1 ^max 1 


Pf^^" (GeV) 


e,/x 


2.4 


5 


jets 


3 


40 


6-jets 


2.5 


40 



Table 3: Acceptance of the different final states in the simulated detector. 

• j"*): m jets and n leptons (electrons or muons) are required in the acceptance 
region with the isolation corresponding to ^ i. 

• CnZ'- n Z boson(s) are reconstructed from lepton kinematics. Lepton candidates 
fulfilling the acceptance cuts must have the same flavour, opposite charges, and a 
Pt > 10 GeV to reduce the amount of leptons from B meson semi-leptonic decays. 
A Z boson is then reconstructed if the di-lepton invariant mass lays in a 10 GeV 
mass window around the Z mass. 

• Cb- at least one of the jet passing the acceptance cuts is 6-tagged. 

• CAih,h)- The two leptons li and I2 have different flavours and opposite charges, 
belong to the same AR < 1.2 cone, and have an invariant mass smaller than 25 GeV. 

In the following, we present simple strategies that can lead to promising Signal-over- 
Background (S/B) ratios. Our purpose is to illustrate the new possibilities that open up 
in the M2HDM and motivate more detailed studies. To this aim detailed information on 
the efficiencies and the visible cross sections are given. The possibility for additional, more 
sophisticated, selection methods is also briefly addressed. 

4.3 hh^H^ ^ ZA^ 

In a type II 2IIDM (e.g., the MSSM scalar sector), the cross section of hh is enhanced 

as tan (3 increases. This process has been shown to offer a promising discovery channel at 
the LHC when the Higgs boson decays into a t^t~ pair (e.g., see Ref. [^] and references 
therein). For the mass spectrum defined by the benchmark point BP2 of the M2IIDM, a 
particularly interesting decay mode is — > Z^ (see Fig. |^. 




Figure 9: Feynman diagram for the hb ZA^ process. 
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The production cross section of bb — > H is about 15 pb when normahzed to the NNLO 



value reported in |77]. Thanks to this large value, a wide range of interesting final state sig- 
natures can be considered. The cases where the ^4" boson provides two 6-jets or r-jets are 
characterised by large rates, yet call for advanced identification and reconstruction tech- 
niques that need a dedicated experimental effort. For the sake of illustration we consider 
here the decay t~^t~ e^ji^ and Z leptonic decays. Requiring different flavors in 

the A^ decay considerably suppresses backgrounds with Z/7* — > e^e~ or fJ,'^fi~. 

The selected final state contains the three same fiavour leptons e"'"e~e^^^ + 
fi'^ fi^ fi^e"^ + ^rp . Such a multi-lepton final state is extremely clean and does not suffer 
from jet reconstruction uncertainties. The relevant backgrounds are Z(Z/^* tt), tiZ 
and W^W~Z with the decay of W and Z bosons into e, fi, t and r ^ e, + ^j-. The cross 
sections for the signal and background processes after all decays are given in Tab. ^. The 



Process 


Decay (MC) 


a X BR{ih) 


ZA 


{Z ^ l+l-){A ^TT^ e±/xT^T) 


4.2 


Z{Zh) 


{Z M-){Z/Y -^TT ^ i+i'-$t) 


10 


ttz 


{t l+b$TW l~b$T)){Z M-) 


3.5 


w+w-z 




0.4 



Table 4: Cross sections of signal ZA^ and backgrounds processes taking into account the leptonic 
final state considered in the analysis. The notation I includes only e and /i, whereas / also contains 
r decaying into e or fi. All cross sections correspond to the final states in the second column. 



selection proceeds as follows. The acceptance cut A{lf) is applied, followed by the Ciz cut. 
The same fiavour opposite sign leptons paired whose mass is the closest to the Z mass is 
retained. The two remaining leptons h and I2 are then assumed to come from the light and 
boosted pseudoscalar, and forced to satisfy CAihih)- The relative and total efficiencies, 
as well as the visible cross sections for all processes listed in Tab. |^ are reported in Tab. |5|. 





ZA 


z{z/r) 


ttz 


Ailf) 


(%) 


51 


18 


42 


Cz 


(%) 


74 


63 


60 


Ca 


(%) 


85 


3.6 


3.3 


^tot 


(%) 


32 


0.39 


0.84 




ifb) 


1.5 


0.039 


0.029 



Table 5: Relative efficiencies of the considered cuts, together with the total efficiencies after all cuts 
and corresponding visible cross sections for signal and background processes. The WWZ process is 
omitted since its visible cross section is four orders of magnitude smaller than that of the signal. 



As can be clearly seen, this channel is very promising: the S/B is high enough so 
that an excess over the SM could be identified after a few inverse femtobarns of integrated 
luminosity. Note that the {Z — > + jets background has also been considered due 
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to its very large cross section {0{nh)), the possibility for jets to produce fake electrons, 
and the possible presence of leptons from heavy meson decays. An inclusive sample of 
10^ events was generated using the matching procedure |^8[ and no event has passed the 
isolation cuts. This background is therefore neglected. However a more detailed study 
should be performed with a more realistic detector simulation and event reconstruction. 

Besides a pure counting experiment a more exclusive study can also be attempted. The 
mass of the two neutral resonances could be measured with an accuracy depending mostly 
on the reconstruction quality. In the signal, the main source of missing transverse 
energy originates from the r's. If the direction of the is required to lay between the 
transverse position of the two leptons li and I2 and the condition > 50 GeV imposed, 
then a proper reconstruction of the invariant mass my^o can be achieved (see Fig. ^). A 
substantial improvement of the S/B ratio is also gained. Finally, the mass can be 
estimated from the and Z boson 4-vectors. 




in-' ' .... I .... I .... I .... I .... I .. I I I I I I I I I I I I I I I I I I I I I I 

^" 0.5 1 1.5 2 2.5 ) 20 40 60 80 100 

A R Reconstructed Mass oflsl' 

Figure 10: Left: the number of events in function of the distance AR between the two leptons not 
assigned to the Z, after applying the A{ll) and Cz cuts. Right: the number of events in function 
of the reconstructed mass of using the leptons 4-vectors and the missing Et after applying the 
additional Ca cut (which require AR < 1.2). Both figures correspond to an integrated luminosity 
of 30 fb^^. The markers show a hypothetical event excess for this luminosity. 

4.4 g{b/b) {t/t)H^ W-{b/h)W+A^ 

The previous analysis shows that, for the benchmark point BP2, a discovery could be made 
after a few inverse femtobarns of integrated luminosity together with the identification of 
two neutral Higgs bosons and ^4'^. However, in order to fully determine the structure 
of an extended scalar sector, it is also crucial to observe a charged Higgs boson. In the 
M2HDM, we expect it to be nearly degenerate in mass with as a consequence of the 
twisted custodial symmetry. 

The associated production of a charged Higgs with a top quark, g{b/b) {t/t)H^ 
(see Fig. |l^), is in general considered as a challenging channel at the LHC. The discovery 
potential strongly depends on tan /?, the mass of the charged Higgs boson and the considered 
decay mode. However, at variance with models such as the MSSM, the M2HDM offers the 
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possibility for the to decay into M^^yl'^. Its observation would therefore be a very strong 
evidence that the scalar sector originates from the M2HDM. We consider this possibility in 
the benchmark point BP2. To normalize the expected signal, we use the NLO prediction 
for the charged Higgs production cross section from Ref. [^], i.e., 465 fb. 



b 




Figure 11: Representative Feynman diagram for the gb tH^ W^W bA'^ process. 



As in the previous analysis, we focus exclusively on the decay —>■ t~^t~ where the 
r+r^ pair decays into e^/j,^. Despite the fact the total signal cross section is reduced 
by almost two orders of magnitude compared to the — > bb case, a strong reduction of 
the background is foreseen if one of the W bosons decays leptonically. The considered 
final state is therefore l^jjbe^fi^ + ^x- If the light quark pair comes from the W boson 
produced in the charged Higgs decay, the resulting jets tend to be collinear due to the large 
boost. As a consequence, they might not be resolved but merged into a single "large" jet 
(noted J) by the reconstruction algorithms. We include both possibilities. 

The relevant backgrounds are W'^W-W^jj, tt{Z/-f*), W{Z/j*)jj, Z{Z/-f*)jj, and 
tW{Z/j*), with j standing for all light and b quarks. The cross sections for the signal 
and the considered background processes, as well as the corresponding final states, are 
summarized in Tab. |6|. For W^/Z + rr + jj and W^W~^W~jj, the jets are initially 
produced with a minimal Pt of 10 GeV, a maximal pseudo-rapidity of 5, and an angular 
separation of AR{jj) > 0.1 for the firsts and AR{jj) > 0.2 for the latter. The details of 
the decay modes and the corresponding rates are shown in Tab. |6[ 

In order to increase the S/B ratio, the acceptance cut A{lf,j'^) and the Cf, cut are 
applied. These are followed by a veto on the presence of a Z boson Cz- The two closest 
leptons {li, I2) with opposite charges and different flavours are assume to come from the 
light and boosted pseudoscalar Higgs boson ^4*^, and therefore are required to satisfy the 
Cyl(^i)^2) cut. The relevance of this last cut is illustrated in Fig. |l^, where the left-hand 
side plot shows the di-lepton invariant mass after applying all cuts except 6*^(^1,^2)- The 
relative and total efficiencies, as well as the visible cross sections for all processes listed in 
Tab. 1^, are reported in Tab. 0. 

A rather low visible signal cross section confirms that this channel is also very chal- 
lenging with the unusual — > VF^^'' decays. However, the S/B ratio of order 0(1) 
leaves some hope that a charged Higgs could still be discovered after a large integrated 
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Process 


Decay (MC) 


ax BR (fb) 


tt{Zh*) 

w{z/r)jj 

Z{Zh)jj 

mw{z/r) 

W+W-W^jj 


{{t/t) ^ l^b$T){H^ ^ {W^ - jj){A ^TT^ e^p.^$T)){b) 
{[i/tj ^ JJ0){H [W L ^TjyA^TT^e fi^^T)j{b) 
{t incl.){t incl.){Z/Y ^ I'+l"' $t) 

iw^ ^ i^$^){z/t - M"'$T)n 

{t ind.){l^ 1'+1"-$t) 

{w^ i^$t){w^ i^$t){w- i'$T)n 


0.75 

4.5 
48 
10 
0.6 
13 



Table 6: Cross sections of the signal gb tH^ — > W^bW^A^ and the relevant background 
processes, taking into account Icptonic and jet final states considered in the analysis. The notation 
I means that the three flavour of leptons are taken into account, and the tau leptons decay into 
e or fi. On the contrary, 7 means that the tau leptons decay inclusively {Z/{Z/'-f*) case). All the 
quoted cross sections correspond to the final states in the second column. 
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Figure 12: Left: invariant mass distribution of the two closest leptons (with different charges and 
fiavours) for both signal and background events. Right: two dimensional distribution of events 
after A{li,j^), Cb, Cz and Ca cuts, as a function of the acoplanarity between the sum of colhnear 
leptons and the hardest non 6-tagged jet acoplanarity aco(2/, J), and the third lepton aco(2/, /) . 



luminosity (~ 300 fb~^). In any case it should be kept in mind that the benchmark point 
BP2 is not the most optimistic scenario: a lighter associated with a larger tan /3 would 
sizably increase the production cross section. 

In addition, more exclusive discriminant variables could be used to exploit further the 
characteristics of the typical topology. As an example, let us consider the fact that the 
heaviest particle in the process is the charged Higgs boson with at least twice the mass of 
the top quark. As a result, it is typically produced with a small transverse momentum, 
giving acoplanar and bosons with large boost. This acoplanarity {A<p between 
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(%) 


■35 


lb 


5.5 


6.0 


14 


21 


Cbi%) 


(%) 


39 


48 


6 


6.3 


39 


49 


Cz{%) 


(%) 


98 


98 


91 


60 


98 


95 


Ca{%) 


(%) 


84 


15 


19 


12 


11 


6.5 


^tot 


(%) 


11 


1.1 


0.056 


0.017 


0.61 


0.64 




Uh) 


0.083 


0.051 


0.027 


0.0017 


0.0037 


0.083 



Table 7: Relative efficiencies of the various cuts together with the total efficiencies after all cuts 
and corresponding visible cross sections for signal and background processes. 



considered final states) can be estimated from the two coUinear lepton e^fi^ together with 
the decay products of the W originating from the charged Higgs. This decay product 
is either the third lepton, or the "large" jet J if the W boson from the charged Higgs 
decays hadronically. Since the two topologies are a priori not known, the two acoplanarity 
definitions (resp. aco(2Z,/) and aco(2/, J)) are built for each event. The distribution of 



signal events with respect to these two variables is illustrated on the right plot of Fig. 12. As 
the distribution of background events is much more uniform in this plane, an enhancement 
of the S/B ratio of around 10% can be achieved if a cut aco(2/, /) + aco(2/, J) > 3 is applied. 

4.5 gg^h^ ^ H^H^ ZA^ZJ^ 

In the context of SM Higgs searches at the LHC, it has been shown that the discovery of 
at high mass could be achieved after a few inverse femtobarns, notably with the observation 
of /i*' — > ZZ or W~^W~ decays [p0|, Bl|. A deviation from the expected visible cross section 



could reveal the presence of additional decays of the Higgs boson such as those predicted 
by the M2HDM. In this context, it is interesting to study the process gg ^ ^ H^H^ 
(see Fig. ^ since it benefits from a relatively large gluon fusion production cross section 
(around 10 pb at m/jO of 400 GeV, when using the NNLO prediction from |^2|) and a sizable 
H^H^ branching ratio. This possibility is investigated assuming the benchmark point 
BPl of the M2HDM and considering the process gg ^ h° ^ H°H° ZA^ZA^. It also 
provides a clear in the final state signature thanks to the presence of two Z bosons (see 
Fig. 1^), required to decay into e and ^. 

Under this decay hypothesis, the total rate decreases considerably, such that only the 
main A^ bb decay mode can reasonably be retained. At the parton level the signal final 
state is then l~^l~l'^l'~bbbb, with = e^,//^. The 6-quarks produced by the pseudoscalar 
A^ are well separated {AR > 1) but have a low average transverse momentum. In order 
to keep a good efficiency by merging the jets two- by- two for the signal, a large cone size 
radius (~ 1) is used. This leads to a final state with only main two jets in the final state. 
It should be noted, however, that with a heavier A^ the individual detection of all four jets 
could be attempted. 

The main backgrounds to be taken into account are ttZ and ZZjj. The process 
gg ^ ^ ZZ can be neglected as well as W^W~ Z because of their relatively low cross 



-28- 




Figure 13: Feynman diagram for tlie gg ~* H^H'^ ZJ^ZJ^ process. 

section and low probability to provide a 6-tagged jet. For ZZjj the jets are produced with 
a minimal Pt of 20 GeV, a maximal pseudo-rapidity of 5 and a AR{jj) > 0.3. The cross 
sections for the signal and background processes are given in Tab. H] 



Process 


Decay (MC) 


a(fb) 


ZAZA 
ZZjj 

ttz 


(Z ^ l+l-){Z r+l'-)bbbb 

{Z^in-){Z^l'+l'-)jj 

{t l+b$T)){t l~b$T)){Z M-) 


3.2 
16 
3.5 



Table 8: List of processes considered in the analysis of the ZJ^ZAl' channel. The notation I 
means that only electron and muons are considered. If the notation I is used instead, all flavours 
are included and the taus are decayed in e or /x. 

In order to increase the ratio of the SjB ratio, the acceptance cut A(/q, j'^) and the 
cut are applied. The efficiency of the is assumed to be the same as for single b quark 
induced jets. We then apply C^z cut, where the invariant mass of the two pairs of same- 
flavour leptons are the closest to the actual Z mass. The relative and total efficiencies for 
all processes listed in Tab. ^, as well as their visible cross sections, are reported in Tab. |9|. 
The visible cross section around 0.3 fb and a. S/B ratio close to 3 suggest that, using only 
the simple discriminant variables described here above, the evidence of such a signal could 
be reached with a total integrated luminosity smaller than 30 fb~^. 





ZAZA 


ZZjj 


ttz 




(%) 


27 


11 


18 




(%) 


50 


7.9 


54 


C2Z 


(%) 


72 


75 


4.1 


(-tot 


(%) 


9.6 


0.63 


0.4 




ifb) 


0.32 


0.1 


0.014 



Table 9: Relative efficiencies (in percent) for each cut presented in the text. Combined efficiencies 
and resulting visible cross sections after all cuts are also shown. 
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Figure 14: Left: the difference AMh between m^^o and the mean of m^fo after A{lQ,j'^), Cb and 
C2Z cuts. Right: invariant masses of the pseudoscalar bosons (two entries per event). The 
markers show a hypothetical event excess for this luminosity. 



After application of the reviewed "standard" cuts, the S/B ratio could be further 
improved by applying a more sophisticated cut, taking advantage of the fact that invariant 
mass of the can, in principle, be fully reconstructed. First, the two masses are 
reconstructed, each from one Z and one jet (which for the signal, has to be understood as 
the single jet induced by the boson decay), such that the difference between the two 
possible mass combinations must be minimal and smaller than 100 GeV. Furthermore, 
the mass difference between the candidate directly reconstructed from the sum of all 
4-vectors of the four leptons and the two jets, and the mean of candidate masses, must 



be smaller than 400 GeV. The relevance of this cut is illustrated in Fig. 14 where the 
left-hand side plot shows the distribution of signal and background events as a function of 
the difference between m/jO and the mean of mji^o, after applying all other cuts. The S/B 
ratio could be easily increased up to 5, to the price of a lower signal visible cross section 
(by roughly 30%). 

The distribution of the reconstructed mass for the events passing all cuts is also 



shown in Fig. 14. It illustrates the possibility to measure this parameter with the simple 
algorithm described above. The resolution could certainly be improved as well as the 
signal significance if a proper jet reconstruction with an optimal cone size and the tracker 
information were taken into account. The determination of the invariant mass of the 
pseudoscalar could also be attempted. 

Note however that the conclusions of this section has to be interpreted with some 
caution since 6-tagging efficiency and jet kinematics have to be re-evaluated when two soft 
and/or collinear 6-induced jets are merged into one jet. This question can only be precisely 
addressed with a proper full simulation of detector effects. 



5. Conclusion 



We have examined a minimal extension of the SM scalar sector based on natural symme- 
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tries. Once a suitable mechanism is provided to generate the required triplet mass splitting, 
the proposed model (which we have dubbed Minimal 2HDM) may fulfill existing theoreti- 
cal and experimental constraints, while giving rise to unexpected phenomenology at high 
energy colliders. 

Our starting point, the generic 2HDM, is a simple yet rich framework which may arise 

in the context of various BSM theories. Genuine SM properties give grounds for specific 
global symmetries, like the custodial and CP ones, which may or may not be satisfied by 
new interactions. Their implementation in the context of the 2HDM sheds new light on 
their possible interplay and allows us to naturally introduce a "twisted" scenario. In this 
new scenario, the successful p 1 phenomenological relation is ensured by the degeneracy 
of a pair of charged scalars {H^) with a scalar and not with a pseudoscalar {^) at 

variance with the usual case in the 2HDM. 

Surprisingly enough, this seemly mild difference opens a window to novel and, to a 
large extent, unexplored phenomenology of the Higgs sector. Due to its vanishing coupling 
to a pair of gauge boson, the pseudoscalar state might have escaped the LEP II searches. 
The mass of charged Higgs boson, on the other hand, is already strongly constrained by 
indirect measurements. In the twisted scenario we propose, the large mass splitting between 
(H^,H^) and is naturally reconciled with tight electroweak precision constraints. A 
SM Higgs boson mass larger than ~ 200 GcV may also be accommodated through small 
breaking of the custodial symmetry, driven for example by loop corrections. 

We have illustrated how an inverted mass spectrum m^o ^ nifjo < m^± ^ ruf^o, 
as compared for example to the typical MSSM scenario, with a relatively light A'^ (e.g., 
nij^o ~ 30 GeV) could meet all the theoretical and present indirect and direct constraints. 
It could even appear as favoured by some precision measurement like the S parameter or 
the muon anomalous magnetic moment. This particular mass spectrum leads to interesting 
new signal opportunities, mainly related to the opening up of new decay channels such as 
hP H^H^,H+H-, iJ± ^ W'^A^ or ZA^. To our knowledge, these signals have 

been only partially covered by analyses in the context of more restricted BSM models like 
the NMSSM or the CP violating MSSM. 

A first analysis, including a fast simulation of detector effects, has been performed for 
three characteristic channels at the LHC. In our scenarios the SM-like Higgs boson is at 
400 GeV and could be easily found in the ZZ and WW decay channels after a few inverse 
femtobarns of integrated luminosity. At that point, the possibility of an extended scalar 
sector and in particular of a M2HDM with an inverted hierarchy should be explored. Our 
results suggest that the neutral Higgs bosons and A'^ could be discovered and their 
mass measured already in the low luminosity phase of the LHC running. As a particularly 
illuminating channel we have proposed to look for production in bb fusion and its 
subsequent decay ZA^ with A^ t'^t~ . In this case a discovery could be achieved 

after only a few femtobarns of integrated luminosity. As an alternative or in addition 
to the above process, we have also studied the gg ^ ^ H^H^ — > ZAPZAP channel 
with yl^ bb, in a Type I scenario. Wc find that despite its rather low cross section the 
signal could be disentangled from the (extremely low) SM backgrounds and could provide 
a rather unique chance to test the inverted mass spectrum hypothesis. Imagining that 
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three new scalars were observed, the search for the last component of the scalar doublet, 
i.e., the charged Higgs, would become the first priority. Such a discovery, however, will 
probably be as challenging as in other more conventional scenarios. As a promising and 
unconventional channel we have considered the associated production of a charged Higgs 
with a single top quark, with the decays — > l^^A^ and — > r+r^. Our analysis 
shows that the rather small signal rate over a large three- VF's background would probably 
need a very large integrated luminosity before being significant. 

The perspectives for extending the present work can be categorised along two main 
directions. On the theoretical side, the twisted version of the custodial and CP symmetries 
can probably be generalised to other sectors of the theory, like the Yukawa sector, or to 
more complex scalar sector extensions, i.e., models with more than two Higgs doublets 
or involving higher representations. Applications to more "ambitious" BSM scenarios like 
Technicolor-inspired models could also be envisaged. 

On the experimental side, more detailed studies would certainly be welcome. Our 
preliminary results reveal some unusual and specific challenges characteristic of proposed 
scenarios. Among others, the identification of highly collinear lepton or 6-jet pairs is 
particularly interesting and appears as a recurrent requirement for several processes. For 
the latter case various innovative solutions to these problems may exist or have already 
been proven to exist (e.g. in [^3|) and some of them might be applied to the study of the 
M2HDM in a near future. 
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